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ABSTRACT 

The recurrent nova T Pyx underwent its sixth historical outburst in 2011, and became the subject of an intensive 
multi- wavelength observational campaign. We analyze data from the Swift and Suzaku satellites to produce a detailed 
X-ray light curve augmented by epochs of spectral information. X-ray observations yield mostly non-detections 
in the first four months of outburst, but both a super-soft and hard X-ray component rise rapidly after Day 115. 
The super-soft X-ray component, attributable to the photosphere of the nuclear-burning white dwarf, is relatively 
cool (~45 eV) and implies that the white dwarf in T Pyx is significantly below the Chandrasekhar mass (~1 Mq). 
The late turn-on time of the super-soft component yields a large nova ejecta mass Mq), consistent with 

estimates at other wavelengths. The hard X-ray component is well fit by a ~1 keV thermal plasma, and is attributed 
to shocks internal to the 2011 nova ejecta. The presence of a strong oxygen line in this thermal plasma on Day 194 
requires a signihcantly super-solar abundance of oxygen and implies that the ejecta are polluted by white dwarf 
material. The X-ray light curve can be explained by a dual-phase ejection, with a significant delay between the first 
and second ejection phases, and the second ejection finally released two months after outburst. A delayed ejection 
is consistent with optical and radio observations of T Pyx, but the physical mechanism producing such a delay 
remains a mystery. 
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1. INTRODUCTION 

The five thermonuclear explosions of T Pyxidis observed in 
1890, 1902, 1920, 1944, and 1966 earned the system its place 
as the prototypical recurrent nova, but have also highlighted 
our poor understanding of many aspects of binary evolution 
and nova theory. The community has waited anxiously for the 
next outburst of T Pyx in order to study this peculiar system 
with modern multi-wavelength capabilities, and T Pyx finally 
obliged by entering its sixth recorded outburst in 201 1 April. 
High-quality panchromatic observations are now revealing a 
host of new surprises for this system. 

A nova is a transient event marking a thermonuclear runaway 
on the surface of an accreting white dwarf. The white dwarf 
accretes hydrogen-rich material from a companion star, and this 
accreted material settles down into a thin degenerate layer on 
the surface of the white dwarf. The pressure and temperature in 
this layer increase until explosive nuclear burning begins, and 
the bulk of the accreted envelope is expelled from the white 
dwarf at hundreds to thousands of km s“ ^ . 

Novae are expected to recur on an accreting white dwarf 
with a timescale primarily determined by the white dwarf 
mass and accretion rate (e.g., Yaron et al. 2005; Wolf et al. 
2013). Predicted recurrence timescales vary widely (~1-10^ yr; 
Yaron et al. 2005), and those novae repeating on historical 


timescales have been dubbed “recurrent” novae. Theoretically, 
we expect recurrent novae to occur in binaries where massive 
white dwarfs accrete at high rates, because more massive 
white dwarfs have higher surface gravities, meaning that the 
critical conditions are reached for smaller accreted envelopes, 
and higher accretion rate systems accrue this trigger mass in 
less time. 

Even before 201 1, T Pyx flew in the face of our expectations 
for recurrent novae. The evolution of its optical light curve is 
slow, showing a several months-long plateau around maximum 
light and a relatively slow decline from this maximum (Schaefer 
2010). T Pyx has a short orbital period (1.83 hr; Uthas et al. 
2010), solidly below the cataclysmic-variable (CV) period gap. 
According to the theory of CV evolution, such short-period 
systems should have, on average, very low accretion rates (e.g., 
Knigge 2011), but observations in quiescence — and the short 
nova recurrence time — imply that T Pyx has an accretion rate 
orders of magnitude higher than these expectations (Gilmozzi 
& Selvelli 2007; Selvelli et al. 2008). In addition, measured 
binary parameters imply that the white dwarf in T Pyx may be 
significantly less massive than the Chandrasekhar mass (Uthas 
et al. 2010), in contrast with common assumptions for recurrent 
novae. 

The high accretion rate in T Pyx does not appear to be 
sustainable, as it is exceeds expectations by several orders 
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of magnitude for mass transfer rates driven by gravitational 
radiation (the commonly accepted mass transfer mechanism at 
such short orbital periods). Knigge et al. (2000) and Schaefer 
et al. (2010) have hypothesized that mass transfer in T Pyx is 
in a short-term elevated state, perhaps incited by a powerful 
nova outburst which occurred during the 1 800s (before regular 
records were kept on T Pyx). Before this postulated event, T Pyx 
may have been a typical CV below the period gap, with a very 
low accretion rate and long intervals between novae. However, 
after the hypothesized explosion, the hot white dwarf irradiated 
the companion star and induced an unusually high mass transfer 
rate. Perhaps this irradiation power is slowly dwindling and the 
accretion rate is gradually declining, explaining the increasing 
intervals between nova events observed for T Pyx throughout 
the last century. Observational tests of this hypothesis have 
reached divergent conclusions as to whether there is evidence 
for a secular decline in T Pyx’s accretion rate (e.g., Schaefer 
et al. 2013; Godon et al. 2014). 

Regardless of T Pyx’s history, it is clear that T Pyx pro- 
vides an opportunity to test an unusual corner of nova param- 
eter space. Compared with other novae, we have a thorough 
understanding of the binary system’s parameters and the accre- 
tion rate (Selvelli et al. 2008; Uthas et al. 2010). With multi- 
wavelength data collected from the 2011 outburst, we can mea- 
sure key properties of the nova event, like ejected mass, and 
compare them with predictions from nova models. With this 
goal in mind, campaigns have been carried out across the en- 
tire electromagnetic spectrum, providing an exquisitely detailed 
picture of the 2011 nova outburst of T Pyx (Chesneau et al. 
2011; Kuulkers et al. 2011a, 2011b; Shore et al. 2011, 2013; 
Evans et al. 2012; Imamura & Tanabe 2012; Nelson et al. 2014; 
Ederoclite 2013; Williams 2013; Schaefer et al. 2013; Patterson 
et al. 2013; Sokoloski et al. 2013; Tofflemire et al. 2013; De 
Gennaro Aquino et al. 2014; Godon et al. 2014; Surlna et al. 
2014). 

At every wavelength studied so far, the 2011 outburst of T 
Pyx shows surprising features when compared to expectations 
for recurrent novae. As in previous outbursts, the optical light 
curve shows a sort of plateau for three months (Figure 1), 
implying that the optical photosphere is roughly constant in size 
for ~90 days after thermonuclear runaway (Shore et al. 2013). 
Optical spectroscopy shows that T Pyx is an unusual “hyper- 
hybrid” nova, switching from He/N class to Feu class around 
Day 10, and then back to He/N on Day 65 (Williams 2012, 201 3; 
Ederoclite 2013; Surina et al. 2014). The radio light curve and 
optical measurements of the change in the binary period after 
the nova event imply a large ejected mass (~10“"^-10“^ Mq\ 
Nelson et al. 2014; henceforth. Paper I; Patterson et al. 2013), 
rather than the ~10“®-10“’ Mq expected for recurrent novae. 
Radio light curves also show a late and steep rise (Figure 1; 
Paper 1), implying that either the ejecta in T Pyx were very 
cold (<200 K) during the first ~50 days of the outburst, or the 
bulk of the nova ejecta stalled at an ~AU-scale radius until it 
was finally expelled ~50 days after the thermonuclear runaway 
(such a delay might also explain the long plateau in the optical 
light curve; Shore et al. 2013). 

In this work, we focus on X-ray observations of the 2011 
outburst of T Pyx obtained with Swift and Suzaku, and compare 
our results with inferences from other wavelengths. X-ray 
emission from novae can be split into two broad classes, which 
may, but need not, exist contemporaneously (Krautter 2008). 
The first class is super-soft X-ray emission, characterized by 
effective temperatures between 10^ and 10^ K, and luminosities 
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Figure 1. Overview of the 2011 outburst of T Pyx at optical, X-ray, and radio 
wavelengths. Top: V-band optical data from the AAVSO. Middle: 0.3-10 keV 
X-ray light curve obtained with Swift /'KKY. Bottom: VLA radio light curve at 
2.5,7, and 37 GHz (Paper I). 

(A color version of this figure is available in the online journal.) 

in the range 10^®-10^^ erg s“' . High-resolution spectra obtained 
with the grating instruments onboard Chandra and XMM- 
Newton have confirmed that this emission originates near the 
white dwarf photosphere (Nelson et al. 2008; Rauch et al. 2010; 
Ness et al. 2011; Orio 2012). In some novae (e.g., V2491 Cyg, 
RS Oph) the soft X-ray flux is continuum emission that most 
likely originates at the white dwarf photosphere. In other cases, 
the soft X-rays are associated with strong lines of H and He- 
like carbon, nitrogen and oxygen that likely indicate scattered 
photospheric emission (see Ness et al. 2013 and references 
therein). Super-soft X-ray emission only becomes visible at 
later stages of the nova outburst, once the ejecta have become 
optically thin to the radiation from the hot, still burning white 
dwarf surface layers; therefore, the emergence time of the super- 
soft source can be used as a diagnostic for ejecta mass (Henze 
et al. 2011; Schwarz et al. 201 1). 

Observations of novae in the ~1-10 keV energy range also 
reveal harder X-ray emission on timescales of days to years 
after outburst (e.g., Mukai et al. 2008). During the novae in RS 
Oph and V407 Cyg, both of which have red giant secondaries, 
hard X-ray emission was detected at early times and attributed 
to the interaction of the nova ejecta with the dense wind of 
the companion (Sokoloski et al. 2006; Nelson et al. 2012). 
In systems with less evolved donors (and hence lower density 
circumbinary environments), internal shocks within the ejecta 
have been proposed as the origin for hard X-ray emission (see 
Mukai et al. 2008 and references therein). In these cases, hard 
X-ray emission can tell us about the structure of the ejecta 
as the nova outburst progresses (O’Brien et al. 1994). In the 
last three years, novae have been identified as a new class of 
GeV gamma-ray transient by Fermi/LAT, indicating that the 
shock interactions in novae are capable of accelerating particles 
to relativistic speeds (Abdo et al. 2010; Hill et al. 2013; see 
also Tatischeff & Hernanz 2007). Observations of novae in 
the 1-10 keV range are required to characterize these shocks 
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and fully understand the gamma-ray production mechanism. We 
note that T Pyx was not detected with Fermi (C. Cheung 2013, 
private communication), making it a useful comparison case for 
the study of why some novae produce detectable gamma-rays 
while others do not. 

In this paper, we discuss both modes of X-ray emission 
during the 2011 outburst of T Pyx. In Section 2, we discuss the 
X-ray observations and data reduction; Swift monitoring reveals 
the X-ray evolution at high cadence, while our single epoch of 
Suzaku spectroscopy provides high signal to noise on Day 194. 
In Section 3, we describe the Swift X-ray light curve, and 
Section 4 presents our spectral analysis of the Swift and Suzaku 
data. Section 5 analyzes the observed hard X-ray component 
and concludes that it is likely produced by a shock within 
the ejecta, rather than interaction between the nova and pre- 
existing circumbinary material. Section 6 presents a super-soft 
component with a relatively cool temperature and late turn- 
on time (compared to other recurrent novae). In Section 7, we 
discuss how these X-ray results align with optical and radio 
observations, and suggest that all three wavelength regimes 
support a second, massive, delayed ejection in T Pyx. We 
conclude in Section 8. 

Throughout the paper, we take 2011 April 14 (MID = 55665) 
as to or Day 0, the beginning of optical rise and the start of the 
outburst (Waagan et al. 2011; Schaefer et al. 2013). We also 
assume a distance to T Pyx of 4.8 ± 0.5 kpc (Sokoloski et al. 
2013). 

2. OBSERVATIONS AND DATA REDUCTION 
2.7. Swift 

The 2011 outburst of T Pyx was monitored frequently with 
the Swift satellite as part of the ongoing monitoring of novae 
by the Swift Nova CV Group (Kuulkers et al. 2011a, 2011b; 
Osborne et al. 2011a). A series of 221 Swift X-ray Telescope 
(XRT) observations of T Pyx were carried out between 2011 
April 14 (Day 1) and 2012 April 17 (Day 369), resulting in 
one of the most detailed X-ray light curves of a nova obtained 
to date. All observations were made in photon counting mode. 
The total XRT exposure time was ~271 ks, with a median 
snapshot duration of ~1 ks. Swift observed T Pyx daily during 
the first month of the outburst, and then reduced the cadence to 
2-3 times per week for the subsequent three months. Once the 
X-ray emission became bright at the start of 2011 September 
(see Section 3 below), the cadence was once more increased to 
daily, and remained so for the next two months. The observing 
cadence was reduced at late times, ultimately resulting in weekly 
observations after 2012 January, through 2012 April (Days 
262-369). 

We created grade 0-12 event light curves and time-resolved 
hardness ratios using the Swift/XKT products generator devel- 
oped by the Swift group at the University of Leicester. This 
web-based tool can create X-ray light curves, spectra and im- 
ages of any object that has been observed with the XRT, and 
makes use of tools originally developed for automatic reduc- 
tion of gamma-ray burst observations. Details of these products 
are given in Evans et al. (2009). We used an adaptive binning 
strategy to maximize the count rate and hardness ratio informa- 
tion in the light curves. Prior to Day 124, and after Day 200, 
the data were binned to give a minimum of 20 counts per bin. 
In between these days while the source was bright, the data 
were binned by observation (i.e., all snapshots with a single 


http://www.swift.ac.uk/user_objects/ 



Energy (keV) 

Figure 2. X-ray spectrum of T Pyx as observed with Suzaku on Day 194 in 
three different XIS units (colored as black, red, and blue lines). The data are 
shown in the upper panel as points, overlaid with the best-fit model (the sum 
of a blackbody, a thermal bremsstrahlung plasma, and a Gaussian line) in solid 
lines. The parameters of this model can be found in Table 1 . Residuals in the fit 
to the data are shown in units of sigma in the lower panel. The data shown here 
are without the energy offset discussed in Section 4. 1 , to illustrate the offset in 
energy calibration between XISl and XIS0/X1S3. 

(A color version of this figure is available in the online journal.) 

observation ID). We also created spectra for certain time in- 
tervals during the outburst (see Section 4.2), again using the 
XRT products generator. We obtained the appropriate response 
matrix file, in this case swxpc0tol2s6_20010101v013.rmf, from 
the calibration database. The resulting spectra were binned to 
have a minimum of one count per bin, and then modeled in 
XSpec using the Cash statistic, a maximum likelihood-based 
statistic for Poisson data (Cash 1979). 

2.2. Suzaku 

Given the X-ray rise reported by Osborne et al. (2011a), 
we requested and were granted a Director’s Discretionary Time 
observation of T Pyx with the Suzaku observatory. The exposure 
was carried out on 2011 October 25 (Day 194). Here, we focus 
on the data obtained with the X-ray Imaging Spectrometer (XIS) 
in the 0.3-10 keV energy range. All three functioning XIS units 
were operated in the full-window imaging mode, obtaining 
X-ray event data every 8 s over the full 19' by 19' field of 
view. After standard screening,'^ the net exposure time was 
~3 8,600 s for each XIS unit. We extracted the source spectra 
from a 180" region centered on T Pyx using XSelect v2.4b. 
Background spectra were extracted from annular regions also 
centered on the source, with inner radius 240" and outer radius 
420" for both XISl and X1S3, while the outer radius was 390" 
for XISO (the smaller size is to avoid dead regions of the chip). 
We created response files using the xisrmf gen and xisarf gen 
ftools. The most recent version of the contamination model 
(version 20120719) was used to calculate the effective areas. 
The resulting spectra were binned to have a minimum of 30 
counts per bin in order to facilitate the use of statistics in 
determining the best-fit spectral model (Figure 2). 

3. X-RAY LIGHT CURVE OF THE 201 1 
OUTBURST OF T Pyx 

We present the Swift XRT light curve in the middle panel 
of Figure 1, as well as the V-band AAVSO optical light curve 


** http://heasarc.gsfc.nasa.gov/docs/suzaku/processing/criteria_xis.html? 
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Figure 3. (a) Swift/XRT light curve of T Pyx in the 0.3-10 keV energy range (as in Figure 1, this time with a linear time axis). The red arrow marks the time of 
the Suzciku observation, and the dashed vertical lines mark the period for which we analyze the Swift spectra, (b) X-ray hardness ratio (defined as the ratio of the 
0.55-10 keV to 0.3-0.55 keV count rate) as a function of time, (c) X-ray flux in 0.3-10 keV range derived from weekly averaged Swift spectra. Black diamonds show 
the total flux, red squares mark the blackbody component, and blue crosses denote the bremsstrahlung flux, (c) and (e) Temperature and volume emission measure 
evolution of the thermal plasma component, (f) The evolution of the absorbing column density. The range of interstellar absorbing column derived by Shore et al. 
(201 1) is shown as black dashed lines, and the Gilmozzi & Selvelli (2007) measurement is bracketed by red dotted lines, (g) and (h) The evolution of the blackbody 
temperature and luminosity. The best-fit values from the Suzaku spectrum are shown as red diamonds in panels (d)-(h). 

(A color version of this figure is available in the online journal.) 


(upper panel) and a three-frequency radio light curve obtained 
with the Karl G. Jansky Very Large Array (VLA) in the lower 
panel. The Swift XRT light is also plotted in Figure 3(a) with a 
linear time axis. 

Individual X-ray observations during the first 100 days of the 
outburst mainly resulted in non-detections, with the notable 
exception of the detections between Days 14-20. T Pyx is 


detected as a faint hard X-ray source in this third week 
after outburst (Figure 3(b)), but then the X-rays fade out of 
detectability again for the next three months. 

On Day 117, T Pyx was seen at 0.003 counts s“', and by 
Day 125 it had increased in count rate by a factor of 10. T Pyx 
then entered a sun-angle constraint, and was not observed again 
until 2011 September 3 (Day 142). This first visit back to T Pyx 
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Table 1 

Best-fit Model Parameters for Day 194 Suzaku Spectrum^ 



tbnew* (bb+brems+gauss) 

tbnew* (bb+vapec)fl 

tbnew* (bb+vapec)^ 

tbnew* (bb+vnei) 

N(H) (10^* cm-2) 

0.27 ±0.19 

2.0 ± 0.02 

0.61 ± 0.02 

0.5 ± 0.02 

kTBB (eV) 

47 ±2 

35 ± 2 

44 ±2 

46 ±2 

NoriTlBB^ 

U.UUiO_QQQQg 

n ns'?+o.oo2 

U.U8J_QQQJ 

0.0032 ± 0.0020 

U.UUZO_Q QQQ7 

kTbr (keV) 

0.72 ± 0.03 




Normbr‘= (10-3) 

2.5 ± 0.2 




^^Vapec (keV) 


0.70 ± 0.02 

0.73 ± 0.04 


Nomivapec'* (lQ-3) 


0.25 ± 0.03 

7.4 ± 0.5 


z/Zq 


1.0 

<0.003 


O/Oe 


1507^4 

1 6+^-3 


kTvnei (keV) 




0.70 ± 0.02 

Nornivnei"* (lQ-3) 




7 4+1.0 

Z/Zo 




0.0047«o°o“oi 

O/Oo 




0-8^3 

T (10** s cm-3) 




9 q -1-0.8 

^gauss (keV) 

^•O^J_0.003 




Nomigauss'^ (10“3) 

U-oO_oo8 




XIS 1 normalization 

1.0 

1.0 

1.0 

1.0 

XIS 0 normalization 

0.92 

0.92 ± 0.03 

0.92 ± 0.03 

0.92 ± 0.03 

XIS 3 normalization 

0.87 

0.85 ± 0.03 

0.86 ± 0.03 

0.86 ± 0.02 


657.07 

895.58 

665.24 

653.34 

dof 

539 

540 

539 

537 


Notes. 

“ All quoted uncertainties are 90% confidence intervals. 

*’ NonuBB = where L 39 is the bolometric luminosity of the source in units of 10 ^^ erg s“*, and Dio is the distance to the source in units of 10 kpc. 

"10 

Normbr = — f rtgiiidV, where D is the distance to the source in units of cm and and n; are the electron and ion number densities, respectively, in 

units of cm^^ . 

Nonovapec = / itei^idV, where the parameters have the same meaning as in the bremsstrahlung model. 

NonUgauss = Total photons cm“^ s“* in the Gaussian line. 


revealed that the source had increased in count rate by another 
order of magnitude, to around 0.4 cts s“'. The 0.3-10 keV 
count rate remained high, although variable, until Day ~180, 
at which point it began to fade systematically over the fol- 
lowing 200 days. By the time of the last monitoring observa- 
tion with Swift on 2012 April 17, the source count rate was 
~0.07 counts s“^ 

4. X-RAY SPECTRA OF THE 201 1 OUTBURST OF T Pyx 
4.1. A Deep View on Day 194 with Suzaku 

Although the Swift observations were carried out with high 
cadence, in general, the continuous observation durations were 
rather short and resulted in limited counting statistics for 
spectroscopy. The Suzaku observation on Day 194 resulted in 
~22,000 counts between all three XIS units (Figure 2). We use 
this deep observation to characterize the X-ray emission in as 
much detail as possible, and to inform our model for the lower 
count statistic Swift data. 

We modeled the spectra extracted from all three XIS units 
jointly in Xspec version 12.8.0m (Arnaud 1996), and present 
the resulting parameters in Table 1. In all models, a constant 
was included in order to allow for uncertainty in the calibration 
of the relative effective areas of the three CCDs (typically of the 
order of5%-10%). 

To account for absorption by the interstellar medium (ISM) 
we use the tbnew model of Wilms et al. (2000), and assume the 
cross-sections of Verner et al. (1996) and the ISM abundances 
of Wilms et al. (2000). A range of values for the degree 
of interstellar absorption toward T Pyx have been given in 
the literature. Gilmozzi & Selvelli (2007) find a value of 


E(B — V) — 0.25 ± 0.02 from WE ultraviolet spectra of 
T Pyx obtained during quiescence, which corresponds to a value 
of N{H) Rs (1.3 ± 0.1) X 10^' cm“^ (assuming the correlation 
between N{H) and E(B—V) of Predehl & Schmitt 1995). Shore 
et al. (201 1) used diffuse interstellar bands observed during the 
2011 outburst to determine E{B — V) — 0.49 ± 0.17 mag for 
T Pyx, implying N{H) Ri (2.6 ± 0.9) x 10^' cm“^. Finally, 
Godon et al. (2014) determined an intermediate reddening value 
from UV spectroscopy, E{B — V) = 0.35 mag. 

The Suzaku spectrum is clearly complex. A distinct soft 
component is observed below 0.5 keV and is most obvious 
in the XISl data, which has the highest sensitivity at low 
energies. A second, harder continuum component extends out 
to ~6 keV. Finally, there is a rather striking resolved emission 
line at £ Rs 0.65 keV that we identify as O viii Lya. There are 
no obvious emission lines at any other location in the spectrum. 

We begin our exploration of the spectrum with a simple model 
that accounts for the main features seen in the data: the sum of 
a blackbody, a thermal plasma, and a Gaussian emission line 
(model tbnew* (bb+brems+gauss) in Table 1). This model is 
not strictly physical, as the soft component is expected to show 
features typical of a white dwarf atmosphere and the plasma 
is expected to cool primarily through lines, but it does provide 
a useful starting point given the limited energy resolution of 
the spectrum. This simple model provides a reasonable fit to 
the data, with = 1-26 (v = 539). The Gaussian line 
is clearly required to account for the feature at ~0.65 keV. 
However, the best-ht energy of the emission line is found to be 
0.643;!lo QQ ^3 keV, implying that the O viii line is redshifted by 
~5000 km s“'. This velocity is at odds with the results of 
Tofflemire et al. (2013), who reported that the Oviii Fya line 
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is blue-shifted by 400 km s“^ in a Chandra Low Energy 
Grating Spectrograph (LETG) observation taken a few months 
later (Day 210). The observed redshift is also larger than 
velocities derived from optical spectroscopy at any point during 
the outburst. 

A careful examination of the residuals around the emission 
line in all three XIS units reveals a shift to lower energies in 
the XISl data, possibly indicating a problem with the energy 
calibration of the observation. At low energies (<1 keV), 
uncertainties in the energy scale calibration come primarily from 
the determination of the CCD zero-level, and it is possible that 
this was more uncertain for the XISl CCD than the others. 
We investigated the presence of this energy offset in XISl 
relative to the other two CCDs using the gain f it command in 
XSpec. Fitting all three spectra independently with an absorbed 
blackbody model, we find evidence of an offset in the energy 
scale of XISl by approximately — 8 eV relative to XISO and 
X1S3. Including this additional energy offset in the model fit 
decreases both the reduced value of the fit, and the residuals 
around the emission line in the XISl data. We therefore include 
an energy offset of —8 eV for the XISl spectrum in all of 
our fits moving forward. This energy offset is consistent with 
the estimated uncertainty in the energy scale below 1 keV of 
<10 eV,'^ although it may indicate that the error in the energy 
scale determination for this observation was larger than typical 
values. Comparing fits with and without the offset, we find very 
little difference in the values associated with the blackbody and 
bremsstrahlung components, but the energy of the line shifts 
to higher energies with the inclusion of the energy offset. We 
therefore caution that the uncertainties in the line energy are 
larger than the statistical errors from the fit alone, likely of the 
order of 10 eV (or 4500 km s“^ in velocity space). 

The best-fit parameters obtained for the simple black- 
body-tthermal plasma-tGaussian line model, obtained with the 
inclusion of the additional XISl energy offset, are shown in 
Table 1. This model is presented in Figure 2, overplotted on 
the observed spectra. The best-fit N(H) value for this model is 
only (2.7 ±1.9) x 10^° cm“^, much lower than the value ex- 
pected from optical and UV estimates of reddening. The black- 
body temperature, A:7 bb, is found to be 47 ± 2 eV (545,400 ± 
23,000 K, 90% confidence level), with the main constraint com- 
ing from the well-sampled Wien tail between 0.4 and 0.6 keV. 
The normalization of this component implies a blackbody lu- 
minosity of (3.7]|)^j g) X 10^^ erg s“'. The plasma temperature is 
found to be 0.72 ± 0.03 keV, or (8.4 ± 0.4) x 10® K. The nor- 
malization of this component corresponds to a volume emission 
measure of EMy —n^V — (2.2 '^q\) x 10®^ cm“^, or an X-ray 
luminosity of (8.9;!)og) x 10^^ erg s“^ 

We have accounted for the presence of the O viii Lya line in 
this simple model with an additional Gaussian. However it is 
more realistic to attempt to fit the data with a model that produces 
both continuum and line emission. Therefore, we also modeled 
the harder component with the APEC models of collisional 
ionization equilibrium plasmas (Foster et al. 2012). A solar 
abundance model gives a very poor fit to the data — the lines 
of He-like Si and S, as well as the Fe L-shell emission around 
0.9 keV are much too strong in the model compared to the 
data, while the Oviii Lya line is too weak. If we allow the 
abundances of all metals to vary in tandem (while keeping their 
relative abundances at solar values), the model can achieve a 
good fit to the data with no residuals at the Si and S lines if 
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the best-fit metallicity is extremely subsolar (Z Ri 0.01 Zq). 
However, such a model still fails to reproduce the O viii Lya 
line observed in the data, and results in abundances seemingly 
at odds with the report of approximately solar abundances in the 
larger scale T Pyx ejecta (Williams 1982). 

We obtain a better, but still inadequate, fit to the data 
if the abundance of O in the plasma is extremely super- 
solar, while assuming solar abundance for all other elements 
(model tbnew* (bb+vapec)a in Table 1). Leaving both the O 
abundance of the plasma to vary freely and keeping all other 
abundances fixed at solar values relative to hydrogen, we find 
X^/v = 1.7 (v = 540) for an oxygen abundance of 150jl4 times 
solar. The temperature of this plasma is similar to that found for 
the simple bremsstrahlung model (k?br = 0.70 ± 0.02 keV, or 
8.1 X 10® K). The temperature of the blackbody is slightly lower 
(kTsB = 35 eV, ?bb = 406,000 K). The most marked difference 
is in the N(H), which is found to be more in line with the optical 
and UV values, N(H) = (2.0 ± 0.02) x 10^^ cm-^. 

The large oxygen abundance inferred by the strong emission 
line in our spectrum could indicate that the shocked material 
has been significantly polluted by material from the white 
dwarf. Interestingly, Tofflemire et al. (2013) find evidence for 
a strong overabundance (relative to solar) of N in their high 
spectral resolution Chandra /LETG data from Days 210 and 
235; our earlier Suzaku data unfortunately do not have sufficient 
resolution at low energies to explore this finding. Both oxygen 
and nitrogen overabundance are suggestive of the presence of 
dredged-up white dwarf material in the thermal plasma. 

We also tested an alternative possibility: that the X-rays come 
from a primarily oxygen-rich plasma that has very little metal 
content. We fit an absorbed blackbody + vapec model to the 
data, but this time allowed all elemental abundances to vary 
freely (model tbnew* (bb+vapec);, in Table 1). Given the lack 
of obvious spectral features, we kept the relative abundances of 
He, C, N and all metals with higher atomic numbers than O fixed 
at their solar ratios relative to Fe, and allowed the abundance of 
Fe to vary freely. The abundance of oxygen was allowed to vary 
independently. This model results in a significantly better fit to 
the data than the previous model, with x^/v = 1.2 (v = 539). 
The abundances of all metals other than O are essentially zero, 
and the abundance of oxygen relative to hydrogen is ~60% 
enhanced relative to solar values. The temperatures of the two 
components are similar to those found for other models. 

Alternatively, the bright oxygen line could be an indication 
that the plasma is underionized and has not yet come into colli- 
sional equilibrium (i.e., in non-equilibirum ionization; NEI). To 
test this scenario, we use the vnei model in XSpec. The overall 
spectrum is modeled as the sum of an absorbed blackbody plus 
NEI plasma (Model tbnew* (bb+vnei) in Table 1). Keeping 
all elemental abundances at their solar values, we find no fits 
with x^/J* < 1.7 (y = 538). Allowing the oxygen abundance to 
vary freely, but keeping all other elements fixed at solar does not 
improve the fit much (x^/t < L6, v = 537). Allowing all other 
elements to vary in tandem from their solar values, as we did for 
the vapec model above, does result in a better fit {x^ /v — 1.2, 
V — 537), and results in a similar ratio of oxygen to other metals 
(see Table 1 for model parameters). The best-fit ionization age, 
2.9 X 10'^ s cm“^, is close to the equilibrium value for a plasma 
with kT = 0.7 keV (Smith & Hughes 2010), suggesting that 
non-equilibrium ionization physics is not the primary reason 
for the strong oxygen relative to other metals. 

In summary, the key results of our modeling of the Suzaku 
spectrum are (1) all models require both a soft and hard 
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Table 2 

Best-fit Model Parameters for Weekly Averaged Swift X-Ray Spectra^ 


Date Range 

Time since 

N(H) 

^Tbb 

NormBB^ 

kTbr 

Normbr^ 

•^gauss 

normgauss® 

c-stat 

D.O.F. 

(UT) 

Outburst (Days) 

(lO^l cm~2) 

(eV) 


(keV) 

(10-3) 

(keV) 

(10-3) 



2011 Sep 3-9 

142-149 


32+4 

^•^^-0.46 


T 4+0.6 

■^•^-0.5 

U.OOO_QQi4 

1 oq+0.59 
^•■^^-0.33 

268.1 

267 

2011 Sep 10-16 

149-156 

2-27°o'6 

407s 

U.UJJ_QQ4Q 

1 o 1 -h0.20 

^■^^-0.17 

3 1+0.7 
■^•^-0.6 

0.666" 

n 04+0.86 
^•^^-0.60 

220.8 

222 

2011 Sep 17-23 

156-163 

1 ^-fO.5 

457® 

0 01 ^+0-027 

u.uiD_ooo8 

1 17+0.15 
'-0.16 

-T 0+0.6 

■^•■^-0.5 

0.666" 

Q co+0.50 
^•DD_o.3g 

273.7 

243 

2011 Sep 24-30 

163-170 

1 O-I-0.5 

437® 

0.0197°o1'3^ 

1 qq+0.16 

3.47®-| 

0.666" 

Q 7Q+0.64 
^•'^-0.46 

237.7 

232 

201 1 Oct 1-7 

170-177 

2 

^•■^-0,5 

397® 

U.U01_QQ4-7 

1 Q 1+0.25 
^■^^-0.07 

31^6 

0.666" 

1.DU_o7i 

201.78 

215 

2011 Oct 8-14 

177-184 

1 4+0.4 
^•^-0.4 

487®6 

0 007 

U.UU/1_QQQ4 


23.07°o.; 

0.666" 

^•^O_0.35 

229.4 

216 

2011 Oct 16-21 

184-191 

0 7+0-4 

587; 

U.UUiU_QQQQ5 

0 Q2+0.19 

2 q+0.6 

^•^-0.5 

0.666" 

<0.82 

161.0 

189 

201 1 Oct 23-29 

191-197 

1 2+0.9 

537\; 

0 0022'*'^-^^^ 
U.UUZZ_QQQjy 

o.797«o;; 

2 0+1.5 

0.666" 

Q 4Q+0.64 
^•^^_0.38 

129.8 

146 

201 1 Oct 30-Nov 5 

197-204 

1 2+°-^ 
^•^-0.7 

537‘® 

0 001 

U.UUiU_QQQQg 

0.587“o.\®o 

^•°-l.l 

0.666" 

0 39+0.47 
^•■^^-0.26 

122.7 

148 


Notes. 

“ All quoted uncertainties are 90% confidence intervals. 
*’ NonuBB, nonubr, and normgauss as defined in Table 1. 
Parameter was fixed in model fit. 


component to describe the emission, and (2) only a plasma that 
is significantly enhanced in oxygen relative to all other elements 
can reproduce the observed Suzaku spectrum at £ > 0.5 keV. 

4.2. Spectral Fits from Swift 

To examine the time evolution of the distinct components we 
identihed in the Suzaku spectrum, we extracted and modeled 
weekly averaged spectra from the Swift observations. The flux 
level and total number of counts before Day 142 are too low 
to allow useful spectra to be produced. Likewise, we find that 
after Day 206 the weekly spectra have <100 total counts, and 
so no stable model hts are found. We therefore restrict our 
spectroscopic analysis of the Swift data to weekly averaged 
spectra between Day 142 and Day 206 (2011 September 
3-November 6). At each epoch we model the fit with the simple 
empirical tbnew* (bb+brems+gauss) model developed for the 
Suzaku data, i.e., the sum of a blackbody, a Gaussian emission 
line, and thermal bremsstrahlung absorbed by the interstellar 
medium. Visual inspection of the first few weekly spectra shows 
an enhancement of counts at the position of the Oviii line 
found in the Suzaku spectrum, justifying the inclusion of the 
Gaussian line component. We fixed the energy of the Gaussian 
line because of the lower number of counts in the Swift spectra, 
and allowed the line flux to vary freely. Given the uncertainties 
in the Suzaku energy calibration discussed in Section 4.1, we 
fixed the energy of the line to the best-fit line energy found from 
the Swift September 3-9 spectrum (0.666 keV). Initial values 
for all other parameters are taken from the Suzaku best-ht model 
and then allowed to vary freely. 

The resulting best-fit model parameter values and their 90% 
confidence interval uncertainties are shown in Figure 3 and listed 
in Table 2. Where the Swift and Suzaku observations overlap, 
Swift-denwed parameters are consistent with measurements 
from the deeper Suzaku observation (Suzaku parameters from 
the tbnew* (bb+brems+gauss) model shown as red diamonds 
in Figures 3(d)-(h)). The only exception is the N(H) value, 
which is considerably lower for the Suzaku best-fit model than 
the overlapping Swift data. However, we note that the error bars 
shown on this plot are the statistical uncertainties from the model 
fit only — they do not include a larger systematic uncertainty 
related to the choice of model (see Table 1 for range of N(H) 
values found). 


Over the ~60 days where T Pyx is X-ray bright, we see evolu- 
tion in the spectrum and flux that is dominated by changes in the 
super-soft blackbody component (Figure 3(c)). The temperature 
of the thermal plasma component is roughly constant between 
Days 140 and 180, and then declines by a factor of ~2 by Day 
206 (Figure 3(d)). The volume emission measure over the same 
time period does not vary significantly (Figure 3(e)), and so the 
decline in temperature results in a modest drop in flux from the 
hard X-ray component (see blue points in Figure 3(c)). 

The drop in the total flux by a factor of ~2 between Day 
142-206 is mainly attributable to the blackbody component 
(Figure 3(c)). The decline in the super-soft light curve can 
also be seen in Figure 3(h), where we see that the blackbody 
luminosity is lower after Day 175 than before. This decline 
in flux is despite a roughly constant value of the absorbing 
column density (Figure 3(f); consistent with the interstellar 
values estimated by Gilmozzi & Selvelli 2007 and Shore et al. 
2011) and a blackbody temperature evolution which shows no 
signature of decline (Figure 3(g)). 

In all but one of the spectra, the flux of the emission line 
component is non-zero (Table 2). The flux appears to drop 
slightly in the hrst three weekly averaged spectra, and then re- 
brightens to its highest level in the October 1-7 spectrum. The 
line then fades again; we only obtain an upper limit to the line 
flux in the October 16-21 spectrum. The line flux is non-zero in 
the last two epochs, and ~60% smaller compared to September 
measurements. 

We stress that the X-ray properties shown in Figure 3(c)-(h) 
are only obtained during the plateau in the X-ray light curve. 
Days 142-206. The count rates measured before Day 117 and 
after Day ~250 are orders of magnitude lower, and so we can- 
not glean spectral information at this time due to the poor 
signal to noise in these observations. However, we can infer 
something about the spectral properties from the hardness ra- 
tio plot (Figure 3(b)). We choose 0.55 keV as the break en- 
ergy between the soft and hard bands, since our Suzaku spec- 
tral fits imply that all flux above this energy originates in the 
thermal plasma component. At very early times, the spectrum 
appears to be quite hard, particularly at the time of the first 
X-ray detections around Days 14-20. The spectrum then be- 
comes softer as the source gets brighter. Although the error 
bars are quite large, the hardness ratio during the X-ray rise 
is higher than at peak brightness, which suggests that the hard 
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X-ray component rose before the super-soft emission emerged. 
At late times, we observe a slight trend toward higher hard- 
ness ratios, presumably occurring as the super-soft continuum 
fades. 

5. THE ORIGIN OF HARD X-RAY EMISSION IN T Pyx 

The hard X-rays {E > I keV), modeled in the previous section 
as thermal plasma and seen in all Swift /XRT detections, indicate 
the presence of a long-lived shock. Hard X-ray emission in 
novae can originate either in external material that is swept up 
and shocked by the nova ejecta (Sokoloski et al. 2006; Nelson 
et al. 2008, 2012), or in shocks within the ejecta themselves 
(O’Brien et al. 1994; Mukai & Ishida 2001). In this section, 
we consider the possibility that the hard X-ray emission is 
produced by interaction between the nova ejecta and pre-existing 
circumbinary material. We will show that the timing of the hard 
X-ray component rules out plausible models for circumbinary 
material around T Pyx as a significant source of X-rays. We 
therefore propose that an interaction within the 2011 ejecta, 
between two distinct episodes of ejection, explains the hard 
X-rays. 

The temperature associated with the thermal plasma emission 
encodes information about the kinematics of the shock. During 
the first week T Pyx is visible after solar conjunction (Days 
142-148, the hrst Swift observations with good signal to noise; 
Table 2), the best-fit temperature of the hard X-ray emitting 
thermal plasma is Tbr 1.4 x 10^ K. This temperature is 
comparable to those measured in hard X-ray components of 
other novae (Lloyd et al. 1992; Balman et al. 1998; Mukai & 
Ishida 2001), and can be translated into a rough estimate of the 
shock velocity using the Rankine-Hugoniot jump conditions for 
a strong shock: 


Vs 


l6kTbr 
3 fj. niH 


( 1 ) 


Vs _ I Tbi- 

1000km s“^ “ V 1.4 X 10^ K' 


( 2 ) 


Here, /r is the mean molecular weight of the gas and mn 
is the mass of a hydrogen atom. We therefore estimate 
Vs ^ 1000 km s“' for the X-ray emitting shock in T Pyx. 

The maximum velocity of the ejecta is a critical parameter for 
our models of shocked X-ray emission in T Pyx. We proceed 
with our discussion using H/1 measurements of the maximum 
blueshift from Paper 1; these imply a maximum expansion 
velocity of ~1900 km s“^ on Day 2 and ~3000 km s“' on 
Day 69 (see also Surina et al. 2014). 

5.7. Could Interaction with Circumbinary Material 
Produce the Hard X-Ray Component? 

A rich circumbinary medium is clearly present around T 
Pyx, in the form of a spatially resolved, clumpy Ho!H-[N ii] 
remnant (Williams 1982; Sharaetal. 1997; Schaefer et al. 2010). 
Contini & Prialnik (1997) propose that X-ray emission may be 
produced in T Pyx when ejecta from nova outbursts interact 
with this nebula. Balman (2010) and Balman et al. (2012) 
claimed a detection of spatial extension in X-ray images of 
T Pyx in quiescence, as might be expected for such interaction. 
Subsequently, however, this claim was called into question using 
re-analyzed high-resolution Chandra X-ray images (Montez 
et al. 2012). 


Interaction between the 2011 ejecta and the spatially resolved 
nebula is simple to rule out with timescale arguments. The 
Ho!H-[N ii] flux from the nebula peaks at 4"4, or 0.1 pc (Shara 
et al. 1997). All nine of the nebula’s “shells” proposed by Shara 
et al. (1997) have radii > 1", or 6 x 10'® cm. Similarly large 
radii, 2 x 10*^ cm, are expected for the ejecta from the 1966 
outburst, if they expand at 1900 km s“'. Assuming the 2011 
ejecta are traveling at 1900-3000 km s“', an encounter on 
Day 117 would imply that material from the 1966 ejection is 
traveling at just 20 km s“'. In the unlikely case that there are 
ejecta from 1966 traveling at such a low velocity, we would 
expect the shell of 1966 ejecta to be very thick (spanning 
velocities 20-3000 km s“'), and the shock interaction to be 
prolonged over years. Instead, the rise and plateau of the hard 
X-ray light curve is contained to a relatively short period. Days 
~1 17-206 (Figure 3). Reasonable assumptions for the velocity 
and distribution of the 1966 ejecta therefore imply that they are 
much too distant to be encountered by the 2011 ejecta in the 
hrst year. 

In addition to the ejecta from previous novae, there might 
be circumbinary material present at smaller radii, perhaps lost 
through a wind during the quiescent period, 1967-201 1 . Indeed, 
T Pyx seems to have an unusually high accretion rate for a CV, 
so it is not unreasonable to expect such mass loss (Selvelli et al. 
2008; Schaefer et al. 2013; Godon et al. 2014). However, from 
UV spectroscopy during quiescence, there is no evidence of a 
strong or fast outhow from T Pyx (Gilmozzi & Selvelli 2007). 

Given the low inclination of T Pyx (Uthas et al. 2010), perhaps 
the easiest way to “hide” signihcant mass loss from the binary 
during quiescence is to have it leak out of an outer Lagrangian 
point and remain concentrated in the orbital plane; we might 
then expect mass loss to proceed at relatively modest speeds on 
order of the orbital velocity, ~20 km s“' (Uthas et al. 2010). 
The density of this equatorial circumbinary material should be 
highest near the binary, and fall off linearly with distance. This 
conhguration of circumbinary material is in direct conflict with 
the primary late rise of the X-ray light curve. While the faint 
detections between Days 14-20 might be consistent with such 
an interaction, the X-ray non-detections between Days 20-1 17 
imply a cavity surrounding T Pyx, extending between radii 
~4 X lO*"' — 2 X lO'® cm. No simple picture of circumbinary 
material in T Pyx explains a cavity at these radii, surrounded by 
a dense shell. We therefore move on to consider an alternative 
explanation for the hard X-ray component in T Pyx that rises 
four months after outburst: shocks within the 2011 nova ejecta. 

5.2. Hard X-Rays Produced By Shocks Within the Ejecta 

Shocks may be produced within the nova ejecta themselves, if 
the speed of ejection increases with time after the thermonuclear 
runaway (O’Brien et al. 1994). In this case, the shock velocity 
(Vs) is equal to 4/3 the differential velocity between the inner 
and outer ejecta components. In fact, we see increasing outflow 
velocities in T Pyx over the hrst two months after outburst, as 
measured from optical spectroscopy (Paper I, Surina et al. 20 14). 
The difference between the velocity of material ejection on Day 
2 and Day 69, as measured from the H/i prohle ([3000-1900] km 
s“* = 1 100 km s“*), is consistent with the interaction velocity 
estimated from the X-ray temperature (3/4vj ^ 800 km s“*). 

Let us simplify the complex picture painted by optical 
spectroscopy into a cartoon scenario: one shell is expelled 
on Day 0 with a velocity of ~1900 km s“' and a second is 
expelled later on with a velocity of 3000 km s“' . The bright hard 
X-ray component appears in the Swift data on Day 117, and rises 
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Steeply until it reaches maximum on Day 142 and plateaus until 
Day ~206, after which the flux gradually fades. In our model, 
the period of the X-ray rise and plateau corresponds to the time 
it takes for the shock to plow through the first ejection (O’Brien 
et al. 1994); by Day 206, the first ejection would be completely 
shocked. Using these measured velocities and requiring that the 
outermost edges of the first and second ejections catch up with 
one another around Day 206, we And that the second ejection 
was released on Day 75 (in reasonable agreement with the 
increase of H/1 expansion velocity and the rise of radio flux 
density; Paper I, Surina et al. 2014). 

In addition, the hard X-ray component shows a remarkably 
strong Oviii Lya line in the Suzaku spectrum, implying that 
the shocked gas is overabundant in oxygen by a factor of ~ 1-2 
orders of magnitude over solar values, as is commonly observed 
in nova ejecta (e.g., Gehrz et al. 1998). Such an unusually high 
oxygen abundance implies that the shocked material is nova 
ejecta (rather than solar abundance material which may have 
been stripped from the secondary star). 

The neutral hydrogen column density shielding the shocked 
material decreases only a small amount, if at all, between 
Days 142 and 206, and it is not significantly higher than the 
ISM foreground value during this period (Figure 3(f)). Much 
more dramatic declines in N(H) are seen in the hard X-ray 
components of novae like V1974 Cyg and V382 Vel (Balman 
et al. 1998; Mukai & Ishida 2001). T Pyx’s small decrease in 
N{H) implies that the shock is not viewed through a dense 
absorbing screen and is not deeply embedded in the ejecta. 

We conclude that the timing, temperature, and abundances 
of the bright hard X-ray component are all well described by a 
shock within the ejecta of the 2011 outburst. 

6 . A SUPER-SOFT SOURCE EMERGES AROUND DAY 
130 IN THE 201 1 OUTBURST OE T Pyx 

We clearly observe the unveiling of the super-soft X-ray 
emission in the Swift observations of T Pyx. The super-soft 
phase is thought to begin when the nova ejecta have expanded 
enough to become optically thin to X-ray photons emanating 
from the still-burning shell on the surface of the white dwarf 
(e.g., Krautter 2008; Schwarz et al. 2011). We note that this 
super-soft turn-on is an observationally defined phenomenon; 
the hot white dwarf is almost certainly present at earlier times, 
but it is not observed because of the large absorbing column. 

Our analysis of both the Swift and Suzaku spectra of T Pyx 
show clear evidence for a blackbody component with temper- 
ature of 35^5 eV ([4-5] x 10^ K). The average bolometric 
luminosity of this blackbody between Days 143 and 150 is 
(1 - 3^.0 9 ) X 10^^ erg s“*, near the Eddington luminosity and typ- 
ical of a still shell-burning white dwarf (we note that blackbody 
model fits to super-soft spectra are known to underestimate the 
temperature and overestimate the luminosity of the white dwarf 
photosphere, so the uncertainties are likely larger than the sta- 
tistical errors from the model fit; Osborne et al. 2011b). The 
luminosity of the blackbody is variable and appears to decline 
over the period Days 143-206. Similar variations have been 
observed in other novae, such as V2491 Cyg, where the bolo- 
metric luminosity during the super-soft phase did not appear to 
be constant (Page et al. 2010). 

The observed turn-on time of the super-soft source phase 
in novae has been identified as an important diagnostic of the 
ejecta mass, since more massive ejecta will take longer to thin 
out as they expand, eventually reaching a column density where 
soft X-rays produced at the white dwarf photosphere can be 



Figure 4. When the column density of nova ejecta drops sufficiently (to 
~5 X 10^* cm“^; black solid line), the nuclear-burning white dwarf should 
become visible as a super-soft X-ray source. The vertical gray region shows the 
plausible range of tum-on time for the super-soft phase in T Pyx, implying an 
ejecta mass ^ 10 “^ Mq. We overplot column density as a function of time for 
six different sets of ejecta parameters: those assuming Mej = 10~* Mq (blue), 
Mq = 10“^ Mq (orange), and Mej = 10“^ Mq (red). Solid lines in blue, 
orange, and red represent ejecta with Umax = 3000 km s"*, while dashed lines 
are for Umax = 1900 km s“*. 

(A color version of this figure is available in the online journal.) 


transmitted through the ejecta. Schwarz et al. (2011) and Henze 
et al. ( 2011 ) explore the time elapsed between thermonuclear 
runaway and the observed turn-on of the super-soft source for 
large samples of novae. They find that super-soft turn-on time 
is correlated with the Mej/ Umax- Figure 4 shows the expected 
temporal evolution of the column density for shells ejected 
at Day 0 with a range of masses and two possible expansion 
velocities (chosen based on the H/1 line profiles of T Pyx at 
Days 2 and 69, respectively; Paper I). We assume that the 
shells have an r~^ density profile and spherical symmetry that 
arise from a Hubble flow ejection with Umin/fmax = 0.2 (see 
Paper 1 for details, also Seaquist & Bode 2008). We estimate 
that 5 X 10^* cm“^ is the maximum column density for which 
we could detect a super-soft component with temperature kT < 
50 eV (horizontal black line). 

The super-soft source is definitely observed on Day 143, the 
first observation after the solar gap. A careful examination of 
the data obtained on Days 120 and 123 constrain the earliest 
turn-on time to Day 123, when equal numbers of photons are 
detected above and below 0.5 keV, and a hint of an additional 
soft component can be seen in the spectrum. In Figure 4, we 
And that a turn-on time of 123 days corresponds to an ejected of 
mass of >few x 10“^ Mq, assuming that the ejecta have been 
expanding since Day 0. If there is a two-month long delay in the 
ejection of the bulk of the mass, as suggested by the following 
section, the ejecta mass implied by the turn-on time can be 
slightly lower, but still >10“^ Mq. These mass estimates, while 
larger than traditionally assumed for T Pyx because of its short 
recurrence time (e.g., Contini & Prialnik 1997; Schaefer et al. 
2010 ), are compatible with estimates obtained from analysis 
of the 1966 outburst (Selvelli et al. 2008), radio light curves 
(Paper I), and the decrease of the binary period after the 2011 
outburst (Patterson et al. 2013). 

The blackbody temperature associated with the super-soft 
component is much lower than in other recurrent novae, im- 
plying a white dwarf in T Pyx which is significantly below the 
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Chandrasekhar mass. More massive white dwarfs have smaller 
radii and surface areas, and so have higher temperature photo- 
spheres once they become super-soft sources (Sala & Hernanz 
2005; Wolf et al. 2013). The recurrent novae RS Oph and U Sco 
were both observed to have high blackbody temperatures dur- 
ing their super-soft phases (60-90 eV; Osborne et al. 2011b; 
Orio et al. 2013), as was the classical nova V2491 Cyg (Page 
et al. 2010). T Pyx is strikingly different, with a blackbody tem- 
perature of 30-50 eV, just half of what is measured in other 
recurrent novae (see also Tofflemire et al. 2013). Both the lower 
peak temperature and massive ejected shell are compatible with 
a nova outburst on a lower mass white dwarf (~1 Mq\ Wolf 
etal. 2013). 

The turn-off time of the super-soft source is also theoretically 
predicted to inversely scale with the white dwarf mass. The turn- 
off of the super-soft source refers to an abrupt decline in the 
blackbody luminosity, and likely marks the cessation of nuclear 
burning on the white dwarf surface, due to the exhaustion of the 
reservoir of hydrogen remaining after nova outburst. On a more 
massive white dwarf, the nuclear-burning luminosity is higher 
and the reservoir of hydrogen is less massive, implying a shorter 
duration for the super-soft phase (e.g., Sala & Hernanz 2005). 

We find evidence that the super-soft emission in T Pyx began 
to turn off around Day 180 (Section 4.2). This inference from 
the Swift data is in agreement with the report of a very low- 
luminosity blackbody component in the Chandra spectrum 
obtained on Day 210 (Tofflemire et al. 2013). Wolf et al. (2013) 
predict a stable shell-burning super-soft phase for >1000 days 
following a nova on a 1.0 Mq white dwarf; a super-soft duration 
of ~200 days is more consistent with a nova on a 1 . 15 Mq white 
dwarf. The rapid super-soft turn-off time in T Pyx therefore 
implies a somewhat more massive white dwarf, potentially at 
odds with the late observed turn-on time and the low effective 
temperature in the super-soft phase (e.g., Henze et al. 2014). 
However, as we show in the next section, mass loss from T Pyx 
is peculiar and complex; discrepancies in white dwarf mass 
constraints might be resolved if a larger fraction of the hydrogen 
envelope is expelled from T Pyx than theoretically predicted. 

7. X-RAY, OPTICAL, AND RADIO EVIDENCE EOR A 
STALLED EJECTION IN T Pyx 

In Section 5.2, we found that the hard X-rays in T Pyx are 
consistent with an interaction between two discrete shells: one 
which is triggered at Day 0 and expands at 1900 km s“', and 
the other which begins expanding at 3000 km s“^ on Day ~75. 
This picture is also consistent with peculiarities seen at radio 
and optical wavelengths, as we discuss here. 

The shapes of the radio and X-ray light curves are remarkably 
similar (Figure 1), which at first glance is surprising, as they are 
likely emitted by fundamentally different processes and regions. 
Both essentially show non-detections for the first couple months, 
then rise steeply (Figure 1). However, the quasi-coordinated 
rise is not a coincidence if the increases in both the radio and 
X-ray regimes are tied to the expansion of the ejecta. The radio 
flux rises as the optically thick thermal ejecta expand and the 
emitting area increases in size. The expansion of the ejecta also 
leads to a drop in the column density shielding the hot white 
dwarf, revealing its super-soft X-ray emission. The hard X-ray 
emission is produced when fast-moving ejecta plow into slower- 
moving material. Based on arguments in Section 5, this material 
also appears to originate in the nova ejecta. 

A detailed look at the two light curves reveals that the rises 
in the radio and X-ray are not perfectly synchronized. The first 


detection of the main radio rise was made on Day 67, while the 
first significant detection of the main X-ray rise did not take 
place until 1 .5 months later, Day 117. The lag between the hard 
X-ray and radio light curves cannot be explained if they share a 
common source, with the radio emission originating in the same 
shocked plasma as the X-ray. In addition, the radio emission 
during Days 117-164 is far too bright and optically thick, and 
the temperature of the hard X-ray component is far too hot, to 
share a common source in the shocked gas (see the Appendix). 

In Paper I, we showed that the peculiar radio light curve can 
be explained as either (1) a very cold ejection which is suddenly 
heated (requiring an increase in temperature by a factor of 25 
in ~19 days) or (2) a delay in the expulsion of the bulk of 
the mass by ~60 days after the beginning of the optical rise. 
While both scenarios are consistent with the super-soft X-ray 
evolution of T Pyx, cooling and subsequent re-heating of the 
ejecta provides no natural explanation for the presence of a late- 
onset hard X-ray component that appears to be consistent with 
shock emission. Even when the white dwarf’s photoionizing 
radiation can reach the outer ejecta and heat them, the radiation 
is too soft to result in X-rays with £ > 1 keV. In contrast, the 
second scenario (a delay in ejecting the bulk of the mass from 
the binary; Figure 5) can self-consistently explain both the soft 
and hard X-ray components, along with the radio and optical 
evolution. 

Our proposed multistage ejection scenario is illustrated in 
Figure 5. First, a shell constituting a small fraction (<10%) 
of the total envelope mass (shown in blue in Figure 5) is 
promptly expelled at 1900 km s“' around the time of optical rise. 
Meanwhile, the remaining bulk of the accreted envelope (shown 
in red) puffs up into a quasi-hydrostatic envelope surrounding 
the binary for ~60-70 days. This stalled, massive shell, orbiting 
the central binary system as a common envelope-like structure 
with radius <10'"^ cm, explains the bright, roughly constant 
plateau in the optical light curve and the early-time radio non- 
detections (Hachisu & Kato 2004; Paper I). It is interesting to 
note that models of radio light curves imply that a low-mass 
ejection expelled on Day 0 might account for the sole radio 
detection at early times (measured on Day 17 and followed by 
a suite of non-detections on Days 31^8; Paper I), which is also 
coincident with the very early and short-lived detection of hard 
X-rays. 

Optical spectra reveal that T Pyx originally resembled a 
He/N nova but switched to a Fen nova sometime between 
Day 2 and 10 (Ederoclite 2013; Williams 2013; see Williams 
1992 and Williams 2012 for more interpretation of such spectral 
classifications). A wind with a low mass-loss rate but high 
velocity (M„, 10“* Mq yr“\ 2000 km s“') emanating 

from the puffed-up envelope may sustain the observed high- 
velocity wings on optical emission lines during Days ~2-60 
(Surina et al. 2014) while remaining optically thin and below 
the detection limits in the radio on Days 31-48. 

On Day 65, the optical spectrum of T Pyx was observed to 
transition back to a He/N classification (Ederoclite 2013), im- 
plying that the bulk of the ejected mass was now being lost in 
a discrete shell-like event (central panel of Figure 5). Around 
this time, the optical light curve begins to decline steeply and 
monotonically; we interpret this decline as the optical photo- 
sphere finally receding through the ejecta (Figure 1). In addi- 
tion, the Ha line flux begins to systematically decrease (Surina 
et al. 2014), implying a drop in the density of the envelope. 
It is also around Day 65 that we see the radio light curve 
begin to rise, implying that the envelope is now growing in 
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Figure 5. Schematic illustration of mass ejection and interaction in the 201 1 outburst of T Pyx, integrating constraints from optical, radio, and X-ray data. The yellow, 
blue and red shells represent the large-scale remnant in T Pyx, the first, low-mass ejection, and the later, high-mass shell, respectively. The size of the arrows indicates 
the expansion velocity of each shell. The three panels show key components in the system at Day 20 (left). Day 70 (center), and Day 180 (right); see Section 7 for 
more details. 

(A color version of this figure is available in the online journal.) 


size (Paper I). Lastly, in our internal shock model for the hard 
X-ray emission from T Pyx, we calculated that the faster mov- 
ing ejection was likely expelled around Day 75 (Section 5.2; 
not precisely matching transitions at other wavelengths, hut rea- 
sonably close given the uncertainties involved in our interpre- 
tation of the hard X-ray light curve). The agreement in timing 
from optical, radio, and X-ray tracers is remarkable; all im- 
ply that the bulk of the ejecta mass was not expelled from the 
environs of T Pyx until ~2 months after the thermonuclear 
runaway. 

The physical cause of this stalled expansion remains unclear. 
Other novae have shown similar plateaus around optical maxi- 
mum, with the duration of the plateau widely varying between 
sources, from a few to hundreds of days (e.g., Kato & Hachisu 
2011). However, most are not graced with rich multi- wavelength 
data sets like the community obtained for T Pyx in 201 1 , so these 
plateaus have proved difficult to interpret. Friedjung (1992) 
noted that during the plateau in HR Del, the optical spectrum 
showed little evidence of sustained strong mass loss as seen in 
other novae; it was in fact more consistent with an almost station- 
ary atmosphere. The author proposed that the nova outburst in 
this system was initially not strong enough to eject the accreted 
envelope, and instead settled into a quasi-static conhguration 
with little mass loss. No explanation was offered, however, for 
the transition to the subsequent phase of significant mass loss. 
Hachisu & Kato (2004) concurred that plateaus around optical 
maximum could be interpreted as a static phase of the outburst, 
with an inflated but stable white dwarf photosphere formed after 
the onset of the thermonuclear runaway. 

In later work, Kato & Hachisu (2011) proposed that the in- 
teraction of the companion star with the hydrostatic envelope 
could provide an additional source of energy, empowering the 
outburst to transition from a stable to a mass-losing configura- 
tion. Kato & Hachisu (201 1 ) find that this mechanism should be 
most efficient in very short period binaries, like T Pyx, where 
the companion becomes deeply embedded in the nova envelope. 
In addition, Livio et al. (1990) and Lloyd et al. (1997) hnd that 
novae which undergo such a common-envelope-like phase may 
display ring-like morphologies in their ejecta, aligned with the 
orbital plane; Sokoloski et al. (20 13) hnd hints of a ring structure 
in the Ho!H-[N ii] nebula surrounding T Pyx. 


A different explanation for a delayed ejection from T Pyx 
is proposed by Williams (2012); he suggests that a signihcant 
fraction of the mass loss may have its source in the binary 
companion star, rather than the white dwarf. If the nova 
explosion sweeps up or ablates material from the companion 
star, some of the peculiarities of T Pyx, like its hyper-hybrid 
spectral classihcation, could be explained. Detailed modeling 
of the interaction of the nova event with the companion star is 
required to assess the energetics, mass budget, and timescale of 
the interaction and evaluate if mass loss from the companion 
star might help explain the 2011 outburst of T Pyx. 

Even more speculatively, peculiarities like multiple and 
delayed ejections might be explained if the nuclear-burning rate 
on the white dwarf surface fluctuates with time. Novae show an 
astounding diversity of light curves (Strope et al. 2010), some 
of which are difficult to explain with a single thermonuclear 
runaway. Indeed, nuclear-driven X-ray bursts on neutron stars 
sometimes display multiple thermonuclear runaways in the 
course of a single outburst (e.g., in’t Zand et al. 2003). While 
multiple thermonuclear runaways are not theoretically predicted 
in a single nova outburst, they cannot be ruled out by current 
observations. 

Considering all the available data in tandem, we favor a 
delayed ejection as the primary explanation for the multi- 
wavelength evolution of T Pyx, because this scenario can self- 
consistently explain the late and steep rise in the radio and X-ray 
light curves, the slow evolution of the optical light curve, and 
the unusually complex optical spectral evolution of T Pyx. The 
physical mechanism driving this late-time mass ejection remains 
unclear. 

8. CONCLUSIONS 

In our analysis of Swift and Suzaku data covering the 2011 
outburst of the recurrent nova T Pyx, we detect two distinct 
components of X-ray emission: a super-soft component associ- 
ated with the photosphere of the nuclear-burning white dwarf, 
and a hard component associated with shocked thermal gas. 

The super-soft X-ray component becomes detectable between 
Days 123-143, implying an ejection which is surprisingly 
massive for a recurrent nova (>10“^ Mq), but consistent 
with other recent constraints for T Pyx (Selvelli et al. 2008; 
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Patterson etal. 2013, Paper I). In addition, the temperature of the 
super-soft source is relatively cool (30-50 eV), implying that the 
white dwarf in T Pyx is significantly helow the Chandrasekhar 
mass. T Pyx therefore inhabits a poorly explored corner of nova 
parameter space, where its unusually high accretion rate, rather 
than an unusually massive white dwarf, drives it to have a short 
nova recurrence time. 

A hard X-ray component is also detected in all epochs with 
sufficient signal to noise for spectral analysis (Days 142-206). 
Hard X-rays are relatively common in novae, and are postulated 
to originate in shocks internal to the nova ejecta (O’Brien et al. 
1994; Krautter 2008). The hard X-rays in T Pyx, however, are 
accompanied by an unusually strong O viii Lya emission line, 
which can only be fit if the thermal plasma has a highly super- 
solar abundance of oxygen. Such high oxygen abundances in 
the nova ejecta require the dredge-up of significant amounts of 
white dwarf materials during the nova event (contrary to some 
models of recurrent novae; e.g., Starrfield et al. 2012). Similar 
conclusions can be reached from analysis of Chandra LETG 
spectra, which find a factor of ~15 overabundance of nitrogen 
in the ejecta, relative to solar (consistent with enrichment by 
CNO-processed white dwarf material; Tofflemire et al. 2013). 

In addition, we note the appearance of a faint, hard 
X-ray component at early times with an unknown origin. This 
component was only detectable by Swift [XKY in time-averaged 
spectra spanning Days 14-20, and its low signal to noise pre- 
cludes detailed spectral analysis. However, we note a single 
coincident radio detection with the VLA at 33 GHz on Day 
17 (Paper I). The nature of this component remains a mystery, 
although it may be linked to the first low-mass ejection from 
T Pyx on Day 0 (Figure 5). 

As measured from optical emission line profiles, the expan- 
sion velocity of ejecta in T Pyx increases by 50% during the 
first two months of outburst (Paper I, Surina et al. 2014); we 
find that this variation, and subsequent interaction within the 
ejecta, can naturally explain the temperature and light curve of 
the hard X-ray component. In Section 7, we present a cartoon 
picture that can self-consistently explain the soft and hard X-ray 
evolution, the radio light curve, the optical spectral evolution, 
and the optical light curve. A shell of material is expelled on 
Day 0 at 1900 km s“*, and a second episode of mass ejection 
is released on Day ~65 at 3000 km s“' (Figure 5). While the 
first ejection produces clearer signatures at optical wavelengths, 
the second ejection accounts for the bulk of the expelled mass 
(as implied by the small absorbing column in the hard X-ray 
component and the bright but delayed maximum in the radio 
light curve). We do not yet understand the physical mechanism 
that leads to the bulk of the nova mass stalling out and lingering 
around the binary for two months before finally being expelled, 
but T Pyx now joins a significant and growing sample of novae 
which show evidence for complex, multi-phase mass ejection 
(see Fynch et al. 2008; Krauss et al. 2011; Williams 2012 for 
other examples). 


We are grateful to R. Williams, A. Fderoclite, M. Bode, 
R. Smith, and U. Munari for useful conversations. We thank 
the Suzaku mission for the generous allocation of target-of- 
opportunity time to observe T Pyx. We also thank Neil Gehrels 
and the Swift mission team for their support of the target-of- 
opportunity program for this nova. We acknowledge with thanks 
the variable star observations from the AAVSO International 
Database contributed by observers worldwide and used in this 
research. This work made use of the HFASARC archive. 


data supplied by the UK Swift Science Data Centre at the 
University of Feicester, and observations obtained with the 
Suzaku satellite, a collaborative mission between the space 
agencies of Japan (JAXA) and the USA (NASA). The National 
Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. F. Chomiuk is a Jansky Fellow 
of the National Radio Astronomy Observatory. J. Osborne and 
K. Page acknowledge the support of the UK Space Agency. J. L. 
Sokoloski and J. Weston acknowledge support from NSF award 
AST-1211778. 

Facilities: Swift, Suzaku, VFA, AAVSO 
APPENDIX 

PREDICTING RADIO LUMINOSITY OF THE HARD 
X-RAY EMITTING GAS 


Here we investigate if the same thermal gas which emits 
hard X-rays in T Pyx can account for a significant portion 
of the observed radio flux. Take for example Day 156 (2011 
September 17), a date during the optically thick rise of the radio 
light curve when both VLA radio and Swift X-ray observations 
were obtained. The X-ray observations imply a volume emission 
measure of 2.36 x 10^^ cm“^ and a temperature of 1.2 keV 
(1.4 X 10^ K). 

Assuming a spherical shell with radius of 1900 km s“^ x 156 
Days = 2.6 x 10'^ cm and thickness of 10% the radius, this 
EMy corresponds to a number density of n = 3 x 10^ cm“^. 
The path-length emission measure, EM = f n^dl, determines 
the radio optical depth and therefore radio luminosity; here it is 
EM = 2 X 10^ cm“® pc. The optical depth of radio emission is: 




0.08235 


EM \ 
cm~® pc ) 



(Al) 


(Seaquist & Bode 2008). Therefore, the radio optical depth of 
the shocked gas is r„ = 2 x 10“® at 12.6 GHz, much less 
than unity and implying that the radio emission is completely 
optically thin. The radio flux density expected from this shocked 
gas is: 




(1 - e-^^), 
(A2) 


where is the radius of the shock and D is the distance to 
the emitting body, both in cm. From the hard X-ray emission 
measured on Day 156, we therefore predict a 12.6 GHz flux 
density of 0.02 mJy. On that day, we instead measured a 
12.6 GHz flux density from T Pyx of 11 mJy, a factor of 550 
brighter than expected from the X-ray emitting shocked gas. 

Similar results hold between Days 67-164; the shocked 
material that produces the hard X-rays (or upper limits on the 
hard X-rays) is not sufficient to account for the radio emission. 
The bulk of the radio emission must be coming from another 
source. This estimate underlines the well-known result that 
warm 10“* K gas is an efficient emitter at radio wavelengths, 
while the hot >10® K gas which efficiently emits at X-ray 
wavelengths is usually optically thin and relatively faint in the 
radio. 
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